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Abstract 



J> ' A general framework is presented which unifies the treatment of 

wavelet-like, quasidistribution, and tomographic transforms. Explicit 
[ formulas relating the three types of transforms are obtained. 

The case of transforms associated to the symplectic and affine 
groups is treated in some detail. Special emphasis is given to the 
\ properties of the scale-time and scale-frequency tomograms. Tomo- 

| grams are interpreted as a tool to sample the signal space by a family 

Q_i' of curves or as the matrix element of a projector. 

J: 

1 Introduction 

Several types of integral transforms are used for signal processing in 
^ ■ physics, engineering, medicine, etc. In addition to the traditional Fourier 

analysis M, wavelet analysis has been extensively developed in the last two 
decades § || ||. These two types of transforms are linear transforms. In 
addition, the Wigner-Ville quasidistribution @ , a bilinear transform, pro- 
vides optimal energy resolution in the joint time-frequency domain. A joint 
time-frequency description of signals is important, because in many appli- 
cations (biomedical, seismic, radar, etc.) the signals are of finite (sometimes 
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very short) duration. However, the oscillating cross-terms in the Wigner- 
Ville quasidistribution make the interpretation of this transform a difficult 
matter. Even if the average of the cross-terms is very small, their amplitude 
may be greater than the signal terms in time-frequency regions that carry no 
physical information. To profit from the time-frequency energy resolution 
of the bilinear transforms while controlling the cross-terms problem, mod- 
ifications to the Wigner-Ville transform have been proposed. Transforms 
in the Cohen class U make a two-dimensional filtering of the Wigner-Ville 
quasidistribution and the Gabor spectrogram [|lOj] is a truncated version of 
this quasidistribution. 

Recently, a new type of strictly positive bilinear transforms have been 



proposed, namely, the Radon-Wigner transform |]TT| [|T^] or, more gener- 
ally, the noncommutative tomography [|13|] which, in addition to the time- 
frequency domain, also applies to other noncommutative pairs like time- 
scale, frequency-scale, etc. It is this last class of transforms that will be 
called the tomograms in this paper. The tomograms are strictly positive 
probability densities, provide a full characterization of the signal and are 
robust in the presence of noise. 

Developing a general operator scheme, we show how linear transforms like 
Fourier transform or wavelets are related to the quasidistributions and the 
tomograms. Explicit general formulas are derived relating the three types of 
transforms. The time-frequency plane is then briefly discussed because most 
of the material has been treated before Hl3| . Special emphasis is, however, 



given to the transforms associated to the affine group, namely, to the time- 
scale and frequency-scale tomograms. To clarify the physical meaning of the 
tomograms, we also propose an interpretation as a sampling of the signal 
space by families of curves or as the action of a projection operator. 



2 Wavelet-like transforms, quasidistributions, 
and tomograms 

We present a unified general construction of three types of transforms used in 
signal analysis. The first class consists of wavelet-type transforms, the second 
of quasidistributions, and in the third class are the tomographic transforms. 
Quasidistributions are transforms like the Wigner- Ville one § @] or the P- 
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quasidistributions of Glauber and Sudarshan [ II] [ 13] . Husimi-Kano positive 



quasidistributions |16| [0 will be discussed as well. These types of quasidis- 
tributions are unified in the class of s-ordered quasidistributions [|Tj|. 

In quantum mechanics, quasidistributions describe a quantum state in 
terms of phase-space quasiprobability densities. In signal analysis, quasidis- 
tributions describe the structure of analytic signals in the time-frequency 
plane. There also exist quasidistributions characterizing the signal structure 
in the time-scale plane |L9| [^(J [[Hp p2fl . We refer to quasiprobability densi- 



ties because the corresponding functions are not conventional probabilities, 
being either complex or nonpositive. In the case of positive quasiprobabilities 
like the Husimi-Kano function, the two arguments of the function are not 
simultaneously measurable random variables. The corresponding observables 
do not commute and the uncertainty relation prevents the existence of a joint 
distribution function for noncommuting observables. Time i and frequency 
&, time and scale | (uii + iuij , or frequency and scale are common examples 
of such pairs of noncommuting observables. 

The general setting for our construction is as follows. 

Signals f(t) are considered to be vectors | /) belonging to a dense nuclear 
subspace Af of a Hilbert space 7i with dual space Af* (and the canonical 
identification Af C A/"*). {U(a) : a G /} is a family of operators defined on 
A/"*, and a fortriori on Af by the canonical identification Af C Af* . In many 
cases, the family of operators U (a) generates a unitary group. However, this 
is not a necessary condition for the consistency of the formalism, provided 
the completeness conditions discussed below are satisfied. 

In this setting, three types of transforms are defined. Consider a reference 
vector h G Af* chosen in such a way that the linear span of {U(a)h G Af* : a G /} 
is dense in Af*. This means, in particular, that, out of the set {U(a)h}, a 
complete set of vectors can be chosen to serve as a basis. Two of the trans- 
forms considered are: 

W ( f h \a) = (U(a)h\f), (1) 

Q f (a) = (U(a)f\f). (2) 

If U (a) is a unitary operator generated by B ( c?) = a±t + ««2^ and h is a 

(generalized) eigenvector of the time-translation operator, Wt (a) becomes 
a Fourier transform. With the same B (~ct) plus the parity operator, Qf{&) 
would be the Wigner-Ville transform. Similarly, for B (7?) = ct\D + ia2- d \, 
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where D is the dilation operator D = —\ (it-^ + i-^tj, Wj- (a) is a wavelet 
transform and Qf(a) the Bertrand transform. 

We will denote the transforms of the W^-type as wavelet-type transforms 
and those of the Q/-type as quasidistribution transforms. 

In general, if U (a) are unitary operators, by Stone's theorem, there are 
self-adjoint operators B (a) such that 

Wf\a) = (h | e iB ^ | /), (3) 

Qf\a) = (f\e iB ^\f). (4) 

In this case, because B (a) has a real valued spectrum, another transform 
may be defined, namely, 

Mf\X) = (f\5(B(a)-X)\f). (5) 

This is what we call the tomographic transform or tomogram. In contrast to 
the quasiprobabilities, the transform M^\a) is positive and, as we will see 
below, it can be correctly interpreted as a probability distribution. There- 
fore, it benefits from the properties of the bilinear transforms, without being 
plagued by the interpretation ambiguities associated to the quasidistribution 
transforms. 

For a normalized vector | /), 

(/ I /) = 1, (6) 
the tomogram is a normalized function 

J M { f B) (X) dX = 1 (7) 

and therefore, it may be interpreted as a probability distribution for the 
random variable X corresponding to the observable defined by the operator 
B (n). The tomogram is a homogeneous function 

Mf lv) {X) = \p\Mf\pX). (8) 
The three classes of transforms are mutually related 

Mf\x) = ±-jQf B \a)e-* kX dk (9) 
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and 



Qf\a) = I M { f B/p) {X)e ipX dX. 



'f w - J lvl f y. yi )°~ uyi - ( 10 ) 

Wavelet-type transforms, quasidistributions, and tomograms are related by 
the formulas 



Q 



/ 



CXI 



W ( /\a) 



where 



4 /e 



jx 



M {B) (X) - iM)°\X) - M)"\X) + iM) D A \X) 



r(B) 



(B), 



r( B ), 



(11) 

dX, 
(12) 



I fx) =1 h)+ | /); | / s > =| h)- | /); 
\f 2 )=\h)+i\f); \U)=\h)-i\f). 

Another important case concerns operators U(n), which can be repre- 
sented in the form 

U(a) = e ib ^ P h e~ ib{a \ (13) 

Ph being a projector on a reference vector | h). This creates a quasidistribu- 
tion of the Husimi-Kano type 



Hf\a) 



(f I U(a) | /). 



In the following sections, we show how known examples of wavelet-like 
and quasidistribution transforms are described within the framework pre- 
sented above. We will consider quasidistributions such as Wigner-Ville [0 

0- 



We reformulate the stan- 



[H], Bertrand [19|, and Husimi-Kano [|l 
dard wavelet analysis in terms of this general scheme using an operator U(a) 
belonging to the two-dimensional affine subgroup of the symplectic group 
ISp(2,R). Tomographic transform schemes for time-frequency, time-scale, 



and frequency-scale pairs [0] will be studied within the framework of the 
general approach. Inversion formulas are obtained for the tomograms as well 
as the explicit connection of the wavelet transform to the time-scale and 
frequency-scale tomograms. 



3 Time— frequency transforms 

Here we discuss the case where the operator B (if) , ~dt — (jU, v) is 

(a) = lit + uu 
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(with Cj = — i d/dt) or is equal to this one plus a parity operator. The wavelet- 
like transform in this case is just the Fourier transform and we discuss only 
the tomograms and the quasidistributions. 
The tomogram, that is, 



Mf (X, fi,v) = (f\6 (fit + uu-X)\f) 



was shown in 13 to be 



Mf (X,fi,u) 



2ttM 



exp 



ijj,t 2 itX 



2v 



v 



f(t)dt 



(14) 



(15) 



The tomogram ( |l"5"D is normalized if (/ | /) = 1 

J Mf (X,fi,u) dX = 1. 

From the tomogram Mf (X, (i, u), the signal f(t) may be recovered up to a 
phase 

/(t)f(0) = l|Mf ) (l,/i J t)exp[i(x-^)] dXdfji. (16) 
According to the general scheme, the corresponding quasidistribution is 

Qffau) = (f | e iBl8) ^ | f) = J Mf (X^,v)e lX dX, (17) 



or 



dt. 



(18) 



This quasidistribution is called the ambiguity function in the signal process- 
ing literature 

The tomogram ( |I~5"D and this quasidistribution are related by Eqs. @ 
and (0). The tomogram flT5|) is also related to another quasidistribution, 
namely to the Wigner-Ville quasidistribution WV(t, uj) [[7] || by 



Mf (X, fi,u) = J exp [-ik(X -fiu- vt)} WV(t, uj) 



dk du dr 



The Wigner-Ville quasidistribution is given by the formula 



a 



WV(r,a)= f r + - /* t-- e 



it 



- jfiu 



(19) 



(20) 



6 



The unitary operator U(t,Q), which determines the Wigner-Ville quasidis- 
tribution by 

w(r ) n) = </|t^)(T,n)|/>, (2i) 



is 



U (WV \t, Q) = e 2<nt-™) e i 7 r(P+^-l)/2 ) (22) 

the generator being 



^ A 7T (t 2 + LU 2 - l) 

B {wv) (r, Q) = 2tu - 2ttt + — ± '- . (23) 



4 Wavelets and quasidistributions in the affine 
group 

4.1 Wavelets 

The wavelet transform of a signal f(t) is a linear integral transform decom- 
posing the signal into a set of basis functions 

wj A \s,r)= Jf{t)h: tT {t)dt. (24) 

The wavelets h S T (t) are kernel functions generated from a basic wavelet h{r) 
by means of a translation and a rescaling (— oo<r<oo, s>0): 

hs , T (t) = ^h(—). (25) 



,s \ s 
Using the operator 

U {a) (t, s) = exp(zTo)) exp(z'logs.D), (26) 

where 

iJ = - (to; + tut) = cut + - , with u = —i — , 
equation (^) can be represented in the form 

h s , T (t) = U^(r,s)h(t). (27) 
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For normalized h(t) the wavelets h SjT (t) satisfy the normalization condition 

\h 3jT (t)\ 2 dt = 1. (28) 
The basic wavelet (reference vector) may have different forms, for example, 



h(t) 



e iu t e ~t 2 /2 



71 



or 



h(t) = (1 - t 2 ) e-' 2 / 2 



(29) 
(30) 



called the Mexican hat wavelet. 

The inverse of the wavelet transform wi A \s,r) is 



f{t) = N? lwj A \s,r)^=h 



1 , / 1 — t\ dr ds 

s 



with 



One has the property 



J duo, H{uj) = J h(t) e~ iuJt dt. 



J\wj A \s,r) 



2 dr ds 



N h / \f(t)\ 2 dt. 



(31) 



(32) 



(33) 



Let us consider the operator with the parameters jl = (/ii, H2) = (f-i, v) 
being 

H = - — - . (34) 



v = logs, 

We obtain for the unitary operator (^ 



s - 1 



U {A) \t, s) = U {A \n, u) = exp (ifioo + ivD) 
and the operator B(a) becomes 

B[ A) = iilj + uD. 



(35) 



(36) 



According to the general scheme with the operator (|36|) , the wavelet trans- 
form ( f24j ) can be rewritten in the form ([!]) 



w\ A \s,r)^(h\U^{^y)\}). 



(37) 
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The commutation relation for the operators Co and D is 

[lu,D] = -ilu. (38) 

The commutation relations fl38|) define the Lie algebra of the affine group. 
Therefore, the wavelet transform is the nondiagonal matrix element of a 
unitary irreducible representation of the affine group. The parameters fi and 
v are the group parameters and the Hermitian operator B± belongs to the 
Lie algebra of the affine group. 



4.2 Quasidistributions 

The diagonal matrix elements of the irreducible representation determine a 
quasidistribution for the signal f(t) 

Q ( f A \s,r) = (f\e^ + ^\f), (39) 



with and v expressed in terms of shift and scaling parameters by Eq. (|34j ) . 
Defining the action on the vector | /) as 

e iQ>&+»D)/*f(t) = F(ji,u,t), (40) 

the quasidistribution (|39| ) may be rewritten 

Q { f(s,r) = J F*(-v,-v,t)F(v,is,t) dt. (41) 



Using the known kernel (Green function) of the operator fl35|) ||24|| , one ob- 
tains for the quasidistribution ([11]) the following expression in terms of the 
parameters r and s: 

Qf (s, r) = J > {^jpj f(V~s[t + r/2]) dt. (42) 



The wavelet transform (|24|) may also be written in a similar form 



W ( f A \s,T) =Jh* f (^/s [t + r /2]J dt, (43) 



.* ft 



as follows from Eq. (O 
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4.3 Relation of wavelets to time— frequency tomograms 

In view of ([H]) and flT5|), one relates the wavelet transform to the time- 
frequemcy tomogram 



h(t) 



X 



exp |^ (s 2 t 2 + 2str + r 2 ) - l —(st + r)} ds dt dr 



(44) 



The inverse transform reads 



W { f A \ S ,T) 



2vr£> 
x exp < 



J h*{t)Mf\x,yL, st + r 
ifi(st + T 



where 



V 



1 r 

— / M < f S) {X,^,0)e iX dX dfi 
2ix J 



dX dfj, dt, 



1/2 



(45) 



For the Mexican hat wavelet fl30|), with the admissibility condition 



h(t) dt = 0, 



one has the explicit form 



Mf\x,n,v) 
where 

K M (s,r,fi, u,X) 



1 



27rM|iVj 



J wj A \s,r)KM(s,T, [i,v,X) ds dr 



(46) 



'2tz [s - 



l/JS 



3\ -3/ 2 



(sfIT — sX) 2 ijJLS 1 



x exp 



v 2 — ijius 2 v 
(s/ir — sX) 2 ifiT 2 iXr 



- 1 



and 



2{y 2 — ifivs 2 ) 2v v 
N h = I IculV^ du = 1. 
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The inverse transform for the Mexican hat wavelet reads 

W ( f A) (s, r) = J R(s, r, X, //, v)Mf } (X, fi, v) dX dfi du, (47) 

the kernel being 

1 



R(S,T, X,fl, u) 



2nV^/s 



T — V 



cxp 



2 V s J 2 



4.4 Tomograms. Frequency— scale and time— scale 

The tomogram associated to the operator b[ a ^ = jjuj+uD has been computed 



m 



1|. It is 



1 



2w\v\ 
1 



duo — =- exp 

uj>0 \/0J 



-i [ — io logu; 



2ix\v\ 



o;<0 



do; — ;= exp 



—i 



uj log \u>\ ] 



(48) 



/(cj) being the Fourier transform of the signal f(t). 
The tomogram corresponding to the operator 



B (A) = id + vD 



was also computed, namely, 



M 



(At) 



2tx\v\ 



[ dt ^2 exp 
Jt>o \Jt 



2ix\v\ 



t<0 



dt ^ exp 
t\ 



(-*- 


-logt) 






V J 




i 


- - log 


t 




V 





(49) 



The quasidistribution Q\ (//, z/) related to the above tomogram is constructed 



from the affme group. It was discussed in |19| . To compare signal analysis 
based on the time-scale tomograms and based on wavelets, it is useful to write 
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the tomogram M f (s, fi, v) in terms of the wavelet transform W \ (s, r) 



2ir\uN%\ 



+ 



L 



dtdrds , ft — t\ ,. t . 
n- 2 & r ) ex P 

i>0,s>0 V^SS Z 

(it dr ds , ft — t 



i I -t - — log* 



i<0, S >0 ./|t|s S 2 



ft, ( ) W/(s, r) exp 



i(-t- — logltl 

V V 



>(.50) 



The tomographic transform is invertible, that is, the signal may be recovered 
from the tomogram, namely, 



/(*) = 
and 

f(t)f*(t> 



Mf t} (s,/i,l)e 2iS dn ds 



-1/2 



J M ( f At) (s, /x, 1) e 2iS "^ d/i (51) 



M 



f Au) (s,/x, i/) exp ^ i 



47T 2 

x |u;| d// dz/ do; ds. 



t + t' 

s — //a; — vuj V uj[t — t) 



(52) 



5 Meaning of the tomograms 

5.1 Sampling the phase space 

In the time-frequency space, the tomogram 

Mf (X, n,v) = (f\5 (fit + 1/ cj-X)\f) (53) 

is the expectation value of an operator delta-function in the state | /). The 
support of the delta-function in (|53|) is a line in the time-frequency plane 

X = fxt + vlo. (54) 

Therefore, Mf (X, [A, v) is the marginal distribution of the variable X along 
this line in the time-frequency space. The line is rotated and rescaled when 
one changes the parameters fi and v. In this way, the whole time-frequency 
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space is sampled and the tomographic transform contains all information on 
the signal. 

It is clear that, instead of marginals collected along straight lines on the 
time-frequency plane, one may use other curves to sample this space. For 
the tomograms associated to the affme group, one has 

M<*> (S u v) = (f | 5 (/i + - | /) (55) 

and 

M<^> (S 2 , /i, „) = (/ | 5 + - S 2 ) | /). (56) 

The curves in the time-frequency space, defined by 

Si = fit + utu, S 2 = iao + utu, (57) 
are hyperbolas. This becomes clear using the system of coordinates 

q=-L(t-u), p=-L(t + u). (58) 

In the new coordinates, the curves are 



\ (59) 

which are equations for a parametric family of hyperbolas. This means that 
for the tomograms (Si,fi,u) and M^ A ^ (S?., n,v) the probability dis- 

tribution is collected not on straight lines but on hyperbolas. Other general- 
izations are obvious. One might use marginals on ellipses, parabolas, or on 
any other algebraic curves. 



5.2 Operator delta-function as a projector density 

While constructing tomograms, the nonnegative operator 

8{B{a)-X) (60) 

plays an essential role, B(a) being an Hermitian operator with nondegenerate 
continuous spectrum. The random variable X takes values on the spectrum 
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of B(a). Considering a set of generalized eigenstates (in Af*) of B(a), one 
obtains for the kernel 



(Y | 5 (B(a) - X) | Y') = 8{Y' - X) 5{Y - Y') = (Y \ X) (X \ Y'). (61) 
Therefore, we may identify 5 {B(a) — X) with the projector | X)(X \ 

6{B(a)-X)=\X)(X\=P x . (62) 

From this, it follows 

Mf ] = (f \S(B(a)-X) | /) = (/ | X)(X | /) = \{X \ f)\ 2 (63) 

displaying the positivity of the tomogram. This means that there is always 
a basis in M* such that, by projecting on this basis, the tomogram is the 
product of two complex conjugate functions. By a unitary transformation S, 
B(a) may be transformed to 

SB(a)S ] = B\a). (64) 

Then if {| Z}} is the set of (generalized) eigenvectors of B'(a), js^ | Z)\ is 
a set of eigenvectors for B. Therefore, 

Mf\Z) = (f \6(B(a) - Z) | f) = \(Z \ S \ f)\ 2 = (f \ ^ | Z)(Z \ S \ f). 

(65) 

In this case, the operator U(a) in the general scheme described in Sect. 2 
would be 

U(a) = 5 f | Z)(Z | S, (66) 

which, in this case, is not represented as an exponent of an operator. The 
form (0), with the presence of a projector operator P z —\ Z)(Z \, also shows 
the relation of the tomograms to transforms of the Husimi-Kano type. No- 
tice, however, that, for example, the time-frequency Husimi-Kano transform 

Q f (t,u) = \([3\f)\ 2 , P= tJ ^f, (67) 

where | /3) is a coherent state, does not describe a joint probability distribu- 
tion in the time-frequency plane, because time and frequency do not com- 
mute and a joint probability distribution of two noncommuting observables 
cannot exist due to the uncertainty relation. Therefore, the correct way to 
interpret the Husimi-Kano quasidistribution is not as a joint time-frequency 
probability but as a unitarily transformed tomogram. 
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6 Discrete Spectrum 

In the case where the operator B(a) has a discrete spectrum, one uses a 
Kronecker delta-function 8k (B(a) — n) and associates to the operator the 
Fourier integral on the circle. For example, for the number operator 

8 K (a ] a - n) = (** e ^ a ~ n ) dip; n = 0, 1, 2, 3 . . . , (68) 

V ' Z7T JO 

and the matrix elements of this operator in the Fock basis | m) are 

(ml | 8 K (a ] a - nj \ m) = 8 m > m 8 mn ; m, w! = 0, 1, 2, 3 . . . , (69) 
which equal the matrix elements of the projector 



n)(n 



(70) 



This means that, also in the discrete-spectrum case, the Kronecker delta- 
function of the operator B(ot) is reduced to a projector. The tomogram 
associated to the Wigner-Ville function by this method is given by the so- 
called photon-number tomography |25[ |26[ P7fl , i.e., 

w(n,P) = {f | D*(P) | n)(n \ D((3) \ f) 



v/?2 n n! 



J dt e^ t2 ^ +lV ~ 2tlmp H n (t)f{t + V2Re(3 



(71) 



where the complex number j3 is a linear combination of the parameters \x 
and v. 

One may also construct a tomogram using a Dirac delta-function of the 
same operator B(a) 



M(X,n,r)= J f( t ' + kr)r(t"-kr) 

% klT / 2 . 

exp — cot — \ 2t + t 



-ik(X+w/2)_ 



xe -ikQ(t'+t") 



2 , f 2 



+ 



27rsin(A;7r/2) 
it 



sin(/c7r/2) 



(t' + t") 



The inverse of the transform reads 
1 



71 



exp 



-iVL(t + t')+i[X 



71 



dk dt' dt" dt. 

(72) 

dVL dX. 
(73) 
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In this case the Dirac delta-function of B(a)—X is not reduced to a projector 
density. The tomogram (]72f) corresponds to marginals collected from shifted 
circles in the classical phase space. 



7 Remarks and conclusions 

1 - The main result in this work is the formulation of an unified view for some 
linear and nonlinear transforms through the operator formulation developed 
in Sect. 2. The formulation emphasizes the basic unity of these transforms, 
which are related by explicit formulas. Nevertheless, for each particular 
application, one type of transform may be more convenient than the others. 
In particular, when non-ambiguous joint information on noncommutative 
observable planes is desired, tomograms seem to be the most competitive 
type of transforms. 

The formulation applies both to unitary operators or nonunitary ones of 
the form (0). It is also an appropriate framework to construct new trans- 
forms once a particular aspect of the signal is defined and this is expressed 
through the corresponding operator set. 

2 - The operators U(a) may belong to group representations or be op- 
erators of a deformed group. In the cases we have considered in detail, the 
operators belong to the Lie algebra of the group ISp(2,R), which has the 
following six generators: 



When it is unitary, the operator U (jl) may be considered to be an evolu- 
tion operator associated to an Hamiltonian operator formed from the gener- 
ators of the group. Quasidistributions evolve the state |/ > and project the 
evolved state on the initial condition. On the other hand, the tomograms 
collect the probability density on a family of lines in phase space. In the 
cases that were studied, hyperbolas, straight lines and circles were consid- 
ered. Other types of tomograms might be considered, for example, those 
corresponding to parabolas, that is, X = fit + uuj 2 A similar construction 
might be done for other groups and other algebraic structures like quantum 
groups. 



Li — t, L 2 — u, L 3 — 1, 
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